Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Joint multi-ancestry and admixed GWAS reveals the complex genetics behind human cranial vault shape.

The cranial vault in humans is highly variable, clinically relevant, and heritable, yet its genetic architecture remains poorly understood. Here, we conduct a joint multi-ancestry and admixed multivariate genome-wide association study on 3D cranial vault shape extracted from magnetic resonance images of 6772 children from the ABCD study cohort yielding 30 genome-wide significant loci. Follow-up analyses indicate that these loci overlap with genomic risk loci for sagittal craniosynostosis, show elevated activity cranial neural crest cells, are enriched for processes related to skeletal development, and are shared with the face and brain. We present supporting evidence of regional localization for several of the identified genes based on expression patterns in the cranial vault bones of E15.5 mice. Overall, our study provides a comprehensive overview of the genetics underlying normal-range cranial vault shape and its relevance for understanding modern human craniofacial diversity and the etiology of congenital malformations.

The genetic overlap between osteoporosis and craniosynostosis.

Osteoporosis is the most prevalent bone condition in the ageing population. This systemic disease is characterized by microarchitectural deterioration of bone, leading to increased fracture risk. In the past 15 years, genome-wide association studies (GWAS), have pinpointed hundreds of loci associated with bone mineral density (BMD), helping elucidate the underlying molecular mechanisms and genetic architecture of fracture risk. However, the challenge remains in pinpointing causative genes driving GWAS signals as a pivotal step to drawing the translational therapeutic roadmap. Recently, a skull BMD-GWAS uncovered an intriguing intersection with craniosynostosis, a congenital anomaly due to premature suture fusion in the skull. Here, we recapitulate the genetic contribution to both osteoporosis and craniosynostosis, describing the biological underpinnings of this overlap and using zebrafish models to leverage the functional investigation of genes associated with skull development and systemic skeletal homeostasis.

Targeted Sequencing of Candidate Regions Associated with Sagittal and Metopic Nonsyndromic Craniosynostosis.

Craniosynostosis (CS) is a major birth defect in which one or more skull sutures fuse prematurely. We previously performed a genome-wide association study (GWAS) for sagittal non-syndromic CS (sNCS), identifying associations downstream from BMP2 on 20p12.3 and intronic to BBS9 on 7p14.3; analyses of imputed variants in DLG1 on 3q29 were also genome-wide significant. We followed this work with a GWAS for metopic non-syndromic NCS (mNCS), discovering a significant association intronic to BMP7 on 20q13.31. In the current study, we sequenced the associated regions on 3q29, 7p14.3, and 20p12.3, including two candidate genes (BMP2 and BMPER) near some of these regions in 83 sNCS child-parent trios, and sequenced regions on 7p14.3 and 20q13.2-q13.32 in 80 mNCS child-parent trios. These child-parent trios were selected from the original GWAS cohorts if the probands carried at least one copy of the top associated GWAS variant (rs1884302 C allele for sNCS; rs6127972 T allele for mNCS). Many of the variants sequenced in these targeted regions are strongly predicted to be within binding sites for transcription factors involved in craniofacial development or bone morphogenesis. Variants enriched in more than one trio and predicted to be damaging to gene function are prioritized for functional studies.

Collagen has a unique SEC24 preference for efficient export from the endoplasmic reticulum.

SEC24 is mainly involved in cargo sorting during COPII vesicle assembly. There are four SEC24 paralogs (A-D) in vertebrates, which are classified into two subgroups (SEC24A/B and SEC24C/D). Pathological mutations in SEC24D cause osteogenesis imperfecta with craniofacial dysplasia in humans. sec24d mutant fish also recapitulate the phenotypes. Consistent with the skeletal phenotypes, the secretion of collagen was severely defective in mutant fish, emphasizing the importance of SEC24D in collagen secretion. However, SEC24D patient-derived fibroblasts show only a mild secretion phenotype, suggesting tissue-specificity in the secretion process. Using Sec24d KO mice and cultured cells, we show that SEC24A and SEC24B also contribute to endoplasmic reticulum (ER) export of procollagen. In contrast, fibronectin 1 requires either SEC24C or SEC24D for ER export. On the basis of our results, we propose that procollagen interacts with multiple SEC24 paralogs for efficient export from the ER, and that this is the basis for tissue-specific phenotypes resulting from SEC24 paralog deficiency.

Identification of differentially expressed proteins between fused and open sutures in sagittal nonsyndromic craniosynostosis during suture development by quantitative proteomic analysis.

PURPOSE: Nonsyndromic craniosynostosis (NCS), the premature fusion of cranial sutures, results in an abnormal skull shape and is associated with a significant morbidity. Proteomics is a promising tool for disease characterization and biomarker discovery; we aimed to identify biologically relevant differentially expressed proteins for NCS. EXPERIMENTAL DESIGN: Label-based quantitative proteomic profiling using TMT was performed on protein extracted from mesenchymal stem cells, osteoblasts and bone tissue of five open and five fused sutures of sagittal NCS (sNCS) and analyzed using quantitative LC-MS/MS based bottom-up proteomics. Differential protein abundance between open and fused sutures was determined to identify biologically relevant proteins of interest. Proteins were validated in an independent sample set by western blot and immunohistochemistry. RESULTS: We observed 838 differentially expressed proteins between open and fused sutures of sNCS. Decorin, lumican, and asporin were significantly downregulated while COL4A1 and TGFß1|1 were upregulated in fused compared to open sutures. CONCLUSIONS AND CLINICAL RELEVANCE: The majority of significantly differentially expressed proteins between open and fused sutures were observed in the proteomes of osteoblasts suggesting that protein changes contributing to premature sagittal suture fusion occur predominantly at the osteoblast level. Our findings suggest a possible ineffective ECM deposition at the osteoblast cell stage.

Correction: SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism nearBMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure.

Gain-of-function variants and overexpression of RUNX2 in patients with nonsyndromic midline craniosynostosis.

Craniosynostosis (CS), the premature fusion of one or more cranial sutures, is a relatively common congenital anomaly, occurring in 3-5 per 10,000 live births. Nonsyndromic CS (NCS) accounts for up to 80% of all CS cases, yet the genetic factors contributing to the disorder remain largely unknown. The RUNX2 gene, encoding a transcription factor critical for bone and skull development, is a well known CS candidate gene, as copy number variations of this gene locus have been found in patients with syndromic craniosynostosis. In the present study, we aimed to characterize RUNX2 to better understand its role in the genetic etiology and in the molecular mechanisms underlying midline suture ossification in NCS. We report four nonsynonymous variants, one intronic variant and one 18 bp in-frame deletion in RUNX2 not found in our study control population. Significant difference in allele frequency (AF) for the deletion variant RUNX2 p.Ala84-Ala89del (ClinVar 257,095; dbSNP rs11498192) was observed in our sagittal NCS cohort when compared to the general population (P = 1.28 × 10-6), suggesting a possible role in the etiology of NCS. Dual-luciferase assays showed that three of four tested RUNX2 variants conferred a gain-of-function effect on RUNX2, further suggesting their putative pathogenicity in the tested NCS cases. Downregulation of RUNX2 expression was observed in prematurely ossified midline sutures. Metopic sites showed significant downregulation of promoter 1-specific isoforms compared to sagittal sites. Suture-derived mesenchymal stromal cells showed an increased expression of RUNX2 over matched unfused suture derived cells. This demonstrates that RUNX2, and particularly the distal promoter 1-isoform group, are overexpressed in the osteogenic precursors within the pathological suture sites.

A genome-wide association study implicates the BMP7 locus as a risk factor for nonsyndromic metopic craniosynostosis.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10-8): rs781716 (P = 4.71 × 10-9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10-8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10-9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10-8, OR = 0.45; P = 3.31 × 10-8, OR = 0.45; P = 1.09 × 10-8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10-8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10-6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821-55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.

Pegunigalsidase alfa, a novel PEGylated enzyme replacement therapy for Fabry disease, provides sustained plasma concentrations and favorable pharmacodynamics: A 1-year Phase 1/2 clinical trial.

Pegunigalsidase alfa, a novel PEGylated, covalently crosslinked form of α-galactosidase A developed as enzyme replacement therapy (ERT) for Fabry disease (FD), was designed to increase plasma half-life and reduce immunogenicity, thereby enhancing efficacy compared with available products. Symptomatic adults with FD participated in this open-label, 3-month dose-ranging study, followed by a 9-month extension. Three cohorts were enrolled in a stepwise manner, each receiving increased doses of pegunigalsidase alfa: 0.2, 1.0, 2.0 mg/kg, via intravenous infusion every other week. Pharmacokinetic analysis occurred on Day 1 and Months 3, 6, and 12. Kidney biopsies at baseline and Month 6 assessed peritubular capillary globotriaosylceramide (Gb3) content. Renal function, cardiac parameters, and other clinical endpoints were assessed throughout. Treatment-emergent adverse events (AEs) and presence of immunoglobulin G (IgG) antidrug antibodies (ADAs) were assessed. Sixteen patients completed 1 year's treatment. Mean terminal plasma half-life (each cohort) ranged from 53 to 121 hours. All 11 male and 1 of 7 female patients presented with classic FD phenotype, in whom renal peritubular capillary Gb3 inclusions were reduced by 84%. Mean estimated glomerular filtration rate was 111 mL/min/1.73 m2 at baseline, remaining stable throughout treatment. Three patients developed treatment-induced IgG ADAs; following 1 year's treatment, all became ADA-negative. Nearly all treatment-emergent AEs were mild or moderate. One patient withdrew from the study following a serious related AE. Pegunigalsidase alfa may represent an advance in ERT for FD, based on its unique pharmacokinetics and apparent low immunogenicity.

BBS9 gene in nonsyndromic craniosynostosis: Role of the primary cilium in the aberrant ossification of the suture osteogenic niche.

Nonsyndromic craniosynostosis (NCS) is the premature ossification of skull sutures, without associated clinical features. Mutations in several genes account for a small number of NCS patients; thus, the molecular etiopathogenesis of NCS remains largely unclear. Our study aimed at characterizing the molecular signaling implicated in the aberrant ossification of sutures in NCS patients. Comparative gene expression profiling of NCS patient sutures identified a fused suture-specific signature, including 17 genes involved in primary cilium signaling and assembly. Cells from fused sutures displayed a reduced potential to form primary cilia compared to cells from control patent sutures of the same patient. We identified specific upregulated splice variants of the Bardet Biedl syndrome-associated gene 9 (BBS9), which encodes a structural component of the ciliary BBSome complex. BBS9 expression increased during in vitro osteogenic differentiation of suture-derived mesenchymal cells of NCS patients. Also, Bbs9 expression increased during in vivo ossification of rat sutures. BBS9 functional knockdown affected the expression of primary cilia on patient suture cells and their osteogenic potential. Computational modeling of the upregulated protein isoforms (observed in patients) predicted that their binding affinity within the BBSome may be affected, providing a possible explanation for the aberrant suture ossification in NCS.

A variant associated with sagittal nonsyndromic craniosynostosis alters the regulatory function of a non-coding element.

Craniosynostosis presents either as a nonsyndromic congenital anomaly or as a finding in nearly 200 genetic syndromes. Our previous genome-wide association study of sagittal nonsyndromic craniosynostosis identified associations with variants downstream from BMP2 and intronic in BBS9. Because no coding variants in BMP2 were identified, we hypothesized that conserved non-coding regulatory elements may alter BMP2 expression. In order to identify and characterize noncoding regulatory elements near BMP2, two conserved noncoding regions near the associated region on chromosome 20 were tested for regulatory activity with a Renilla luciferase assay. For a 711 base pair noncoding fragment encompassing the most strongly associated variant, rs1884302, the luciferase assay showed that the risk allele (C) of rs1884302 drives higher expression of the reporter than the common allele (T). When this same DNA fragment was tested in zebrafish transgenesis studies, a strikingly different expression pattern of the green fluorescent reporter was observed depending on whether the transgenic fish had the risk (C) or the common (T) allele at rs1884302. The in vitro results suggest that altered BMP2 regulatory function at rs1884302 may contribute to the etiology of sagittal nonsyndromic craniosynostosis. The in vivo results indicate that differences in regulatory activity depend on the presence of a C or T allele at rs1884302.

ISL1 is a major susceptibility gene for classic bladder exstrophy and a regulator of urinary tract development.

Previously genome-wide association methods in patients with classic bladder exstrophy (CBE) found association with ISL1, a master control gene expressed in pericloacal mesenchyme. This study sought to further explore the genetics in a larger set of patients following-up on the most promising genomic regions previously reported. Genotypes of 12 markers obtained from 268 CBE patients of Australian, British, German Italian, Spanish and Swedish origin and 1,354 ethnically matched controls and from 92 CBE case-parent trios from North America were analysed. Only marker rs6874700 at the ISL1 locus showed association (p = 2.22 × 10-08). A meta-analysis of rs6874700 of our previous and present study showed a p value of 9.2 × 10-19. Developmental biology models were used to clarify the location of ISL1 activity in the forming urinary tract. Genetic lineage analysis of Isl1-expressing cells by the lineage tracer mouse model showed Isl1-expressing cells in the urinary tract of mouse embryos at E10.5 and distributed in the bladder at E15.5. Expression of isl1 in zebrafish larvae staged 48 hpf was detected in a small region of the developing pronephros. Our study supports ISL1 as a major susceptibility gene for CBE and as a regulator of urinary tract development.

Genetic advances in craniosynostosis.

Craniosynostosis, the premature ossification of one or more skull sutures, is a clinically and genetically heterogeneous congenital anomaly affecting approximately one in 2,500 live births. In most cases, it occurs as an isolated congenital anomaly, that is, nonsyndromic craniosynostosis (NCS), the genetic, and environmental causes of which remain largely unknown. Recent data suggest that, at least some of the midline NCS cases may be explained by two loci inheritance. In approximately 25-30% of patients, craniosynostosis presents as a feature of a genetic syndrome due to chromosomal defects or mutations in genes within interconnected signaling pathways. The aim of this review is to provide a detailed and comprehensive update on the genetic and environmental factors associated with NCS, integrating the scientific findings achieved during the last decade. Focus on the neurodevelopmental, imaging, and treatment aspects of NCS is also provided.

Muenke syndrome: An international multicenter natural history study.

Muenke syndrome is an autosomal dominant disorder characterized by coronal suture craniosynostosis, hearing loss, developmental delay, carpal, and calcaneal fusions, and behavioral differences. Reduced penetrance and variable expressivity contribute to the wide spectrum of clinical findings. Muenke syndrome constitutes the most common syndromic form of craniosynostosis, with an incidence of 1 in 30,000 births and is defined by the presence of the p.Pro250Arg mutation in FGFR3. Participants were recruited from international craniofacial surgery and genetic clinics. Affected individuals, parents, and their siblings, if available, were enrolled in the study if they had a p.Pro250Arg mutation in FGFR3. One hundred and six patients from 71 families participated in this study. In 51 informative probands, 33 cases (64.7%) were inherited. Eighty-five percent of the participants had craniosynostosis (16 of 103 did not have craniosynostosis), with 47.5% having bilateral and 28.2% with unilateral synostosis. Females and males were similarly affected with bicoronal craniosynostosis, 50% versus 44.4% (P = 0.84), respectively. Clefting was rare (1.1%). Hearing loss was identified in 70.8%, developmental delay in 66.3%, intellectual disability in 35.6%, attention deficit/hyperactivity disorder in 23.7%, and seizures in 20.2%. In patients with complete skeletal surveys (upper and lower extremity x-rays), 75% of individuals were found to have at least a single abnormal radiographical finding in addition to skull findings. This is the largest study of the natural history of Muenke syndrome, adding valuable clinical information to the care of these individuals including behavioral and cognitive impairment data, vision changes, and hearing loss.

Mutations in SEC24D, encoding a component of the COPII machinery, cause a syndromic form of osteogenesis imperfecta.

As a result of a whole-exome sequencing study, we report three mutant alleles in SEC24D, a gene encoding a component of the COPII complex involved in protein export from the ER: the truncating mutation c.613C>T (p.Gln205(∗)) and the missense mutations c.3044C>T (p.Ser1015Phe, located in a cargo-binding pocket) and c.2933A>C (p.Gln978Pro, located in the gelsolin-like domain). Three individuals from two families affected by a similar skeletal phenotype were each compound heterozygous for two of these mutant alleles, with c.3044C>T being embedded in a 14 Mb founder haplotype shared by all three. The affected individuals were a 7-year-old boy with a phenotype most closely resembling Cole-Carpenter syndrome and two fetuses initially suspected to have a severe type of osteogenesis imperfecta. All three displayed a severely disturbed ossification of the skull and multiple fractures with prenatal onset. The 7-year-old boy had short stature and craniofacial malformations including macrocephaly, midface hypoplasia, micrognathia, frontal bossing, and down-slanting palpebral fissures. Electron and immunofluorescence microscopy of skin fibroblasts of this individual revealed that ER export of procollagen was inefficient and that ER tubules were dilated, faithfully reproducing the cellular phenotype of individuals with cranio-lentico-sutural dysplasia (CLSD). CLSD is caused by SEC23A mutations and displays a largely overlapping craniofacial phenotype, but it is not characterized by generalized bone fragility and presented with cataracts in the original family described. The cellular and morphological phenotypes we report are in concordance with the phenotypes described for the Sec24d-deficient fish mutants vbi (medaka) and bulldog (zebrafish).

Sebelipase alfa over 52 weeks reduces serum transaminases, liver volume and improves serum lipids in patients with lysosomal acid lipase deficiency.

BACKGROUND & AIMS: Lysosomal acid lipase deficiency is an autosomal recessive enzyme deficiency resulting in lysosomal accumulation of cholesteryl esters and triglycerides. LAL-CL04, an ongoing extension study, investigates the long-term effects of sebelipase alfa, a recombinant human lysosomal acid lipase. METHODS: Sebelipase alfa (1mg/kg or 3mg/kg) was infused every-other-week to eligible subjects. Safety and tolerability assessments, including liver function, lipid profiles and liver volume assessment, were carried out at regular intervals. RESULTS: 216 infusions were administered to eight adult subjects through week 52 during LAL-CL04. At week 52, mean alanine aminotransferase and aspartate aminotransferase levels were normal with mean change from baseline of -58% and -40%. Mean changes for low-density lipoprotein, total cholesterol, triglyceride and high-density lipoprotein were -60%, -39%, -36%, and +29%, respectively. Mean liver volume by magnetic resonance imaging and hepatic proton density fat fraction decreased (12% and 55%, respectively). Adverse events were mainly mild and unrelated to sebelipase alfa. Infusion-related reactions were uncommon: three events of moderate severity were reported in two subjects; one patient's event was suggestive of a hypersensitivity-like reaction, but additional testing did not confirm this, and the subject has successfully re-started sebelipase alfa. Of samples tested to date, no anti-drug antibodies have been detected. CONCLUSIONS: Long-term dosing with sebelipase alfa in lysosomal acid lipase-deficient patients is well tolerated and produces sustained reductions in transaminases, improvements in serum lipid profile and reduction in the hepatic fat fraction. A randomized, placebo-controlled phase 3 trial in children and adults is underway (ARISE: NCT01757184).

Genome-wide association study and mouse expression data identify a highly conserved 32 kb intergenic region between WNT3 and WNT9b as possible susceptibility locus for isolated classic exstrophy of the bladder.

Bladder exstrophy-epispadias complex (BEEC), the severe end of the urorectal malformation spectrum, has a profound impact on continence as well as sexual and renal functions. It is widely accepted that for the majority of cases the genetic basis appears to be multifactorial. Here, we report the first study which utilizes genome-wide association methods to analyze a cohort comprising patients presenting the most common BEEC form, classic bladder exstrophy (CBE), to identify common variation associated with risk for isolated CBE. We employed discovery and follow-up samples comprising 218 cases/865 controls and 78 trios in total, all of European descent. Our discovery sample identified a marker near SALL1, showing genome-wide significant association with CBE. However, analyses performed on follow-up samples did not add further support to these findings. We were also able to identify an association with CBE across our study samples (discovery: P = 8.88 × 10(-5); follow-up: P = 0.0025; combined: 1.09 × 10(-6)) in a highly conserved 32 kb intergenic region containing regulatory elements between WNT3 and WNT9B. Subsequent analyses in mice revealed expression for both genes in the genital region during stages relevant to the development of CBE in humans. Unfortunately, we were not able to replicate the suggestive signal for WNT3 and WNT9B in a sample that was enriched for non-CBE BEEC cases (P = 0.51). Our suggestive findings support the hypothesis that larger samples are warranted to identify association of common variation with CBE.

Classic bladder exstrophy: Frequent 22q11.21 duplications and definition of a 414 kb phenocritical region.

BACKGROUND: Classic bladder exstrophy (CBE) is the most common form of the bladder exstrophy and epispadias complex. Previously, we and others have identified four patients with a duplication of 22q11.21 among a total of 96 unrelated CBE patients. METHODS: Here, we investigated whether this chromosomal aberration was commonly associated with CBE/bladder exstrophy and epispadias complex in an extended case-control sample. Multiplex ligation-dependent probe amplification and microarray-based analysis were used to identify 22q11.21 duplications in 244 unrelated bladder exstrophy and epispadias complex patients (including 217 CBE patients) and 665 healthy controls. RESULTS: New duplications of variable size were identified in four CBE patients and one control. Pooling of our previous and present data (eight duplications in 313 CBE patients) yielded a combined odds ratio of 31.86 (95% confidence interval, 4.24-1407.97). Array-based sequence capture and high-throughput targeted re-sequencing established that all breakpoints resided within the low-copy repeats 22A to 22D. Comparison of the eight duplications revealed a 414 kb phenocritical region harboring 12 validated RefSeq genes. Characterization of these 12 candidate genes through whole-mount in situ hybridization of mouse embryos at embryonic day 9.5 suggested that CRKL, THAP7, and LZTR1 are CBE candidate genes. CONCLUSION: Our data suggest that duplication of 22q11.21 increases CBE risk and implicate a phenocritical region in disease formation.